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Synchronization, Control and Stability of Fractional 
Order Hyperchaotic Systems in The Context of The 

Generalized Memory 

Vladimirsky E.I., Ismailov B.I. 


Abstract — In the article offered fractional kinetic model of 
networks with generalized memory. On the basis of fractional 
kinetic model network with hyperchaotic systems, embedded in 
a percolation structure, realized task topologically 
synchronization. While tracing control and stability. Criterion - 
“proximity” capture average return time Poincare. Shows a 
visualization of results. 

Index Terms — fractional kinetic model, topological 

synchronization, generalize memory, Poincare return time. 


I. Introduction 

The most control problems of modern philosophy of 
physics are [1, 2]: problem of singularity (space and 
time-mega world); complementarily problem (determinism 
and causality-micro world); the problem of self-organization 
(chaos and order - macro world); problem of the relationship 
of the empirical and the speculative knowledge. 

Accounting and the impact of these problems is 
considerable interest in the context of the analyses and 
synthesis of fractional kinetic Network models 
synchronization, control and stability. 

Naturally that this paradigm is permitted under the 
synergetic approach while remaining within the physics 
Research Area. 

To continue publishing in this article provides a structure 
“Topological synchronization Tracing control stability” is 
permitted in terms of the kinetic model of the network of 
fractional order. 

An important feature of the model is to determine the 
parameters of the percolation at the micro level. Kinetic 
model realized in the space of Poincare fractional order [3]. 

II. Notation and Preliminaries 

A. n-dimensional-fractional-order chaotic system. 

Consider the following n-dimensional-fractional-order 
chaotic system 

D q X = F{X,X 0 ,0 ) (1) 
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where 



X 


2 



eR 


denotes the 


n-dimensional state vector of the original system; X Q - 
represents the system initial state, 

q = (q 1 ,q 2 ,...,qj e R n is a set of fractional order of the 

original system, and 6 = {d x , 0 2 , . . . , 9 D ) r G R D is the 
value of original system parameters. 

Let the fractional-order derivative of the function f(t) in 
the Caputo sense is defined as: 

D“f(t) = J m - q f ( ' n) (t). ( 2 ) 

Here, q is the fractional order, m is an integer that 
satisfies m — \<q<m,f m (t ) is the ordinary m th 
derivative of / , and is the Riemann-Liouville integral 
operator of order |Ll > 0 , defined by 

J " g (i) = 7 ^ - N' g(A h ’ ( 3 ) 

where /"'Q denotes the gamma function. A particularly 
important case in many engineering applications is 
0 < q < 1 . In this situation, Eq. (2) together with Eq. (3) 



1 

r{l-q) 



-x )~ q f(x)dx 


( 4 ) 


The operator Df is often called “ q th-order Caputo 
differential operator” and will be used throughout the paper. 


B. Network model 

In this paper, we consider a dynamical network consisting 
of N identical nodes with nonlinear couplings, in which 
each node is a n -dimensional system obeying the following 
fractional kinetic equations [4] : 

Dlx — Ax + f(x ) (5) 

where x{t ) = (x x {t ), x 2 {t\... X n {t )) r e R n represents the 

state vector, A e R nxn is a constant matrix, and 
f : R n — » R n represents the nonlinear part of the oscillator 

and is assumed to be smooth enough, q = {q l , 
indicates the fractional orders with all q i e (0,l) and 

Dlx=(pi x x x ,Dl 2 x iy . . .,Dl n x n J . We call system (4) a 
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commensurate order system for q } = q 2 = . . . = q n , 

otherwise system (4) is called an incommensurate order 
system. 

The entire network is a system of nN FDEs particular, 
the state equations are: 

D'!x : = Ax, + f(s(t)), i = 1,2, ...,N (6) 

where x G R describes the state of the i -th node and 

is the coupling signal. These unit systems are 
nonlinearly coupled by the function 

s(t) = h(x \,x 2 ,...,x N ) (7) 

where h : R nN —> R n integrates the state of all units 
forming a complex network. It should be noted that the 
integrated function h can be designed in either linear form 


e.g., h{x l ,x 2 ,...,x N ) = — Xj j or nonlinear form 

v N J=1 J J 

(e.g.,h(x l ,x 2 ,...,x N ) = T^ g(x ; ) with smooth 

TV J 

function g : R n — > R n . From the view of synchronization, 

both kinds of function h would achieve network 
synchronization, but with different synchronized state. 
Therefore, for simplicity, we shall later on restrict ourselves 
to the linear case, i.e., to linear integrated function h . For 
such networks, the state equations can be explicitly written as 


D q x i = Ax t + f 


8 


-Vx. 
Ni t J 


i = 1,2,..., N, 


( 8 ) 


where 8 G R is the coupling strength. 

C. Topology of fractional- order space 

Definition 1 . The number is called as a metric order of a 
compact A 

k = lim{— In N A {s)/ In 8) ( 9 ) 

where 8 - the sphere of radius 8 ; N(s) - number of 
spheres in a final sub covering of a set. 

The lower bound of metric orders for all metrics of a 
compact A (called by metric dimension) is equal his 
Lebesgue to dimension. 

However it appeared that the metric order entered in [5], 
coincides with the lower side the fractal dimension of 
Hausdorff-Bezikovich defined in the terms “box-counting”. 

Takes place 

Theorem 1 [5]. For any compact metric space X . 

dimX = inf { lim — ^ £,d ^ — - : d is a metric on X \ 


where 


-log 8 


\U \ : U is a finite open cov ering of 


N £d (x) = i 

X with mesh < £ 

From here (x.d,) - compact fractal metric space with 
dimension d r . 
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Here it is important to note that at the description of 
properties of systems with fractional structure it is impossible 
to use representation of Euclidean geometry. There is a need 
of the analysis of these processes for terms of geometry of 
fractional dimension. 

Remark. In [6] presented results of communication of a 
fractional integrodifferentiation (in Rimana-Liouville or 
Gryunvalda-Letnikov’s terms) with Koch’s curves. 

It is noted that hiunique communication between fractals 
and fractional operators does not exist : fractals can be 
generated and described without use of fractional operations, 
and defined the fractional operator not necessarily generates 
defined (unambiguously with it connected) fractal process or 
fractal variety. 

However use of fractional operations allows generating 
other fractal process (variety) which fractal dimension is 
connected with an indicator of a fractional 
integrodifferentiation a linear ratio on the basis of the set 
fractal process (variety). 

In [6] fractional integrals of Riman-Liouville are 
understood as integrals on space of fractional dimension. 
Thus the indicator of integration is connected with dimension 
of space an unambiguous ratio. 

In this regard consideration of dimension of chaotic 
systems of a fractional order causes interest. So, in [9] was 
noted that dimension of such systems can be defined by the 
sum of fractional exponents E , and 27 < 3 is the most 
effective. 

Let the chaotic fractional system of Lorentz take place [7]: 

d a , x d* , d ! nm 

— x = o(y-x),-y = pX -y-xz,-z=xy-bz (10) 

here cj = 1 0, p = 28, b = 8/3; 0 < a,p,y< 1, r > 1- 

Then fractional dimension of system of the equations (6) 
will have an appearance [9] : 

a + (3 + y = Z (11) 

So, for example, for Lorentz ’s system with fractional 
exponents (X — j3 = y — 0. 99 , effective dimension 
27 = 2.97 . 

This, in the context of fractional dynamics let X - any set 
of nonlinear physical systems, A a - a subset of a set X of 
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systems of a fractional order with memory A a a X . Then a 
triad {x,A a compact fractional metric space with 

dimension X . 

Let’s designate W e (x,d ) . On the basis [8] and 
remarks (x,A a ,X)czW. 

Let’s consider transformation W at an angle of 
communications of average time of return of Poincare (t^ 

with d and “residual” general memory (GM). 

Here g : (x) => d f l : d f => GM, % : (x) => (g,l). 

From here U £ -the generalized compact metric 

space of Poincare with dimension . 

D. Generalized systems with memory 

Let (x, A a , Cl W ; Z the set of all integers [9]: 

= [0,®). K 0 = (-OO.0), Z* = {1,2,...}, and 

Z< 0 = {0,— 1,— 2, — } . 

Definition 2. Let (z , A " , z) d W be GM . 
GM c/?xZ is called a compact generalized memory of 
[4]: 

gm = q >0 (Jq 0 , (i2) 

where: Q >0 = (JIlLiI./) 

7=0 

and: Q <0 = LKk’^-.L^ + 1 ) 

k=l 

for same finite of observed: 

s. <...<s <s 0 =0 = t. <t. <...<t . 

k i 0 0 1 j 

III. Topological synchronization 
U nlike traditional methods of synchronization in [10] 
proposed the concept of topological synchronization of 
coupled chaotic systems. 

Definition 3[10]. Two systems are topologically 
synchronized, if Poincare return times behave a similar way. 

Thus match the dimensions of these two systems - a 
necessary condition for the topological synchronization: well 
indicates the “average similarity” [10]. It is known that the 
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synchronization feature is the preservation of a certain 
frequency: in this case, is the relationship between the 
Poincare return times. 

Invariance of these ratios, are the timing mode. 

A. Synchronization between the two fractional 
hyperchaotic systems 

In the general case, synchronization of chaotic 
fractional-order systems, united has the form: 

d a X , ( . d a Y - x _. M 

—=m—= g( Y) + u (<) c3> 

where a - order of derivative, , a e ( 0 , 1 ] , master’s XeR n 
and slave YeR n systems; f:R n ^R n and 
g : R —> R - vector fields master and slave systems. In 
the general case condition synchronization systems defined 
as: U{t) = {u v . . , i.e. lim t ^ ao ||X - Y|| , where |*| - 

Euclidean norm. 

B. Topological control 

Consider the following general structure of the 
fractional-order nonlinear system under control. 

o D\x{t) = f(x(t)+Bu(t)), (14) 

where u(f)= [w i (f)u 2 (f)...W m (f)] r is m - dimensional 

input vector that will be used and following control structure 
will be considered for state feedback: 

u ( t )= u ,S t ) +u sA t )’ ( 15 > 

where U eii {t) is equivalent control and M (f) is the 

switching control of the system (10). 

With regard to the task the topological control will be 
submitting a number of definitions [11]. 

Definition 4. The system topologically controllable if and 
only if coincides with X on the basis of the criterion metrics 
“proximity” Hausdorff. 

Theorem 2 [12]. Let E and F is compact subset R n , 
£ > 0 . Hausdorf distance H(E, F ) satisfies the relation. 

H(E,F)<£OE^F + £ and Fc£t s, 

where 8 > 0 the allowable threshold. 

Topological control of the formation of a new 
structure made by the algorithm - “a comparison of the 
prototype”. 
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Common challenges in the implementation of 
topological synchronization and topological control is a base 
- definition, theorem and intelligent iterative algorithm. 

Fractional-order iterative learning control scheme is 
given as [9,12]: 

where «,(<) = K (>)~K (>)■ 



Fig. 1. The basic scheme of iterative learning control with 
Y d (t) being the trajectory, U k {t) and Y k (/) the input 
signal. 


Definition 4. Two systems topologically controllable if and 
only if they are synchronized topologically. 

Remark. If you synchronize at the same time a regularity at 
the exit system there is a situation called passive control. 

Otherwise the usual iterative procedure for the organization 
of regular structure on the system output. 

C. Formation of loss memory 

The [9] mentioned that Poincare return time for 
chaotic systems are determined by the well-known fractal 
dimension of the process. The spectrum of the Poincare 
return time on the reversibility is chaotic systems as well 
hence the equivalence between spectrum and distribution 
generalized memory. 

Therefore, memory loss is determined by the 
difference between the global and local fractal dimensions. 

Thus required to resolve the problem in the context of 
the generalized memory. 

D. Poincare recurrence diagram 

Displayed the system on the two-dimensional square 
matrices [N,N] and of formula [13]: 


where N - number of considered (examined) condition X ; 
8 - size of a neighborhood of a point X at the moment i ; 
I'll - norm; 0 (*) - function of Heaviside. 

E. Topological stability of hyperchaotic - order systems 
Determine the stability of the zero solution on the 


system 


dx 

dt 


Q , Q = {co , Qg \j / . 

1 ’ q I n ) n = 0 y q q ' 


Proposition. Let GM gU be structure of 
generalized memory. If there exists a differentiable observed 
V : R — > R + such that the following hold: 

(/) - if trajectory will pass thought the point 0, i.e. 

V(x) < 0, the system is stable with (( r)) and matches 


GM- 

{ii ) - if trajectory will pass below the point 0, the 
system is asymptotically stable with d f ((^)) • 


IV. Main theoretical results 
A. Mathematic model of T -synchronization, tracing control 
and stability with generalized memory (GM). 

Let mathematical model of the network of fractional order 
in the structure “topological synchronization - tracing control 
- stability” of the form 

D q x i = S(x.)+ C(x ) + S(x.)+ r|(x ) + GM(x.) + X(x ) (16) 
where £(x) - synchronization of algorithm i -th systems, 


C(x) - control algorithm i -th systems, 
5(x ) - stability i -th systems, i = 1,2, 


( x \ 

r/ : D q x t v Ajfac — b cp , where frac (x) 

\T J 

the fractional part map of the sawtooth; 

GM(x ) - Generalized systems whit memory: 


is 


m,£. / .,\ 

R. . 1 =0[£. - k -A. ) /, j = i^j.XGU , 

L J 
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K : D q x PL , PL - percolation lattice. 
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Let U e (x,(r)j is tensioned the mathematical of 
the system D q (x ) , i.e. 

D q (x )—»£/. (17) 

So the analyses and synthesis of model (9) components will 
be implemented in context of compact metric Poincare space 

with dimension (x^ . 

In this section consider the structure “topological 
synchronization - tracing control - stability” of the resolution 
on the example of fractional order hyperchaotic systems. 

B. Information processes in the kinetic model 

In the fractional kinetic model of the network an 
important place occupied by the problem of transport 
information. 

In this context arbitrary kinetic structure with 
elements of heterogeneous fractional-order chaotic maps 
embedding on percolation lattice that is carried out the 
operation in order to from structure of relations links. 

Remark. The formal foundation of system embedding 
is the results of modern general algebra on embedding 
algebraic structures in relative simple algebraic structures. 

Definition 5. Embedding kinetic model in 
fl -dimensional lattice is a map in which the vertices (of the 
system) are map in the lattice nodes and edges (connections) 
are on the lattice lines. 

Axiom embedding. Let G = (V , - kinetic 

structure, where V - set of vertices, V = {Ql }, E - set of 
edges E r = \e i }, i = 1,2, . . n . 

Let the set of n -dimensional lattice LxL in E 2 with 
many of nodes, M = j m. j and many links, 

C = { c.}, j = 1,2, 

Let some operator X embedding G = (V,E'^ in 
n -dimensional lattice Z = ( M ,C ) has the form [14]: 

x:G=>Z, (18) 

where 


Then the percolation lattice will be a geometric and 
dynamic implementation of fractional-order kinetic model of 
the network. 

To describe the topology of a set used geometric 
characteristics of the fractional - called index of connectivity 
[15]. 

It is known that the fractal dimension of percolation cluster 
in general is defined as [3]: 

D = d-k ; (20) 

y 

where d - space of a two dimensional flow in E 2 ; /? - 
coefficient for two-dimensional percolation; 

y = 5 / 36 = 0. 1389 ; y = 1/2 - index of correlation. 

However, the real fractal dimension d f is determined by 

the embedding of heterogeneous systems in a percolation 
lattice. 

Remark. In principle, in general, there should be similarity 
of these dimensions, i.e. 

def 

D~d f . (21) 

It is important to note that the cluster will determine it’s 
the throughput, i.e. 

p : |Li — > d f , (22) 

where |Lt - the throughput; p - operator. 

Thus, in the article the task of implementing “topological 
synchronization - tracing control - stability” structure 
fractional percolation lattice of the cases the presence without 
lass of memory and the memory loss. 

C. Algorithm 

Step 1. Let given to the fractional-order hyperchaotic 
system [16]: 

D l =a(y 1 -*i)+>V 
D q =bx x +cx l z l + w,, 

(23) 

Dl =-e x ' y ' ~dz,, 

K =hy. 

Here, a = 10 , b = 40 , c = 2 , d = 2,5 , 


\ V ->M 


G=>Z = < 


E'^-C 


0, otherwise. 


h = —4 , where q is the fractional-order, q = 0.95 . 

(19) Fractional-order hyperchaotic Chen system following 

[17]: 
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=a l (x 2 -x l ) + x 4 , 
d T 2 =yx l -x 1 x 3 + c l x 2 , 
D ‘l; =x 1 x 2 -b l x 3 , 
D*=j V c 3 +d l x 4 , 


where £/, = q 2 = £/ 3 = q A = 0. 95 , the system parameters are 


(«i ,b t ,c\,d t ,y) = (35,3,28,7 ) . 

Step 2. Simulation of the system (23) - (24) according to the 
algorithm [9] . 


Step 3. Let x= a {x n }^ fj is related observable two 

fractional-order hyperchaotic (23) and (24) systems. 

Step 4. In order to achieve the behavior of synchronization 
between two new hyperchaotic systems by using the proposed 
method, suppose the master is X . 

Step 5. The related observable x perturb of sawtooth wave 


rj\D q x i v A frac 


■ + cp 


, where frac (x) is the 


fractional part. 

frac (x) = x — [x] , A is amplitude, T is the 
period of the wave, and cp is its phase. 

y = * + S(*), Y = {y n }^, YeU.Y is slave 


system (fig. 1). 

Step 6. Iterative learning algorithm for topological 
synchronization on schematic “master-slave” with while 
tracking control for hyperchaotic nonlinear fractional-order 
systems shows figures 2 - a,b,c,d,e, f , g . 


Step 7. On the basis theorem and on iterative procedure 
define the effect of the “proximity” capture of the average 
Poincare return time as a criterion for the hyperchaotic 
topologically synchronization systems while tracking control. 



Loss mem Memory 



C 


d 



e f g 

Fig. 2 Visualization of iterative algorithm 
a,b - first iteration: memory allocation and its recurrence plot; 

c,d - next iteration: memory allocation and its recurrence 
plot; e,f - real memory and its recurrence plot; g - Lyapunov 
indicator of chaoticity, lyap.min=-0,2009, lyap.max=2,7177. 


Step 8. Define the Poincare diagram Dx, Dy [13] and 
fractal dimensions: 

d fl =1.2586, d f2 = 1.2034; 

the average Poincare return time as: 

(z\ = 3. 9721428, (r) 2 = 3. 798857. 

Step 9. Effects of the “proximity” define as: 

/ \ prox 

• / (( T )i’( T ) 2 M^( r )i c =( t ) 2 + ^’ 

prox 

(fj 2 Cl + £ , where £ = 0. 1 7 . 

The figure 1, f shows regularity-like of tracking control. 
Step 10. Define Lyapunov stability. 
i ) Let the Lyapunov function is given in the quadratic form 
as: 

vW=|,« ! - 

ii ) We calculate the total derivative of the function 

V(x): 



-140 -40 


60 160 260 360 460 


Loss mem 6 Memory 



a 


b 


dV 

dt 


6V dx _ 
dx dt q 


Q Q< 0 , 


where f2 — > memory function (fig. 3). 
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Step 77. We make embedding fractional kinetic models in 
percolation lattice of rules (18) and (19). On fractional 
percolation lattice is a series connection of heterogeneous 
hyperchaotic fractional mapping (fig.4.). 




Fig.4. Structured fractional kinetic model. 

The size of the lattice 16x16 . Share indignation 
number: 0.763. Number of runs the program: 50. The 
probability of membership node in percolation cluster: 
0.95249. The average cluster size - 4.6076. Probability 
emergence percolation cluster - 0.96. 
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V. Conclusion 

Resolution structure “topological synchronization - control 
- stability” showed method analyses from generalized 
memory. When tracing control process formed the transition 
of “hyperchaos-chaos” confirmed by a positive Lyapunov 
exponent. 
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